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HIGH STIFFNESS SEALS FOR ROTOR CRITICAL SPEED CONTROL 


by David P. Fleming 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135 

ABSTRACT 

An annular seal is analyzed in which the inlet clearance is larger than 
the outlet clearance; the flow path may be either stepped or tapered. Tliis 
design produces radial stiffnesses 1.7 to 14 times that of a constant- 
clearance seal having the same minimum clearance. When sealing high- 
pressure fluids, such a seal can improve rotor stability and can be used to 
shift troublesome critical speeds to a more suitable location. 

NOMENCLATURE 

A,B,E, G coefficients given by eq. (15) 

C clearance of concentric seal 

D seal diameter 

d hydraulic diameter 

e seal eccentricity 

e^ surface roughness 

F radial restoring force 

f , Fanning friction factor 

H seal clearance ratio, Cj^/C 2 

h local seal clearance 

I (l/h^)dz 
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integral defined by eq. (10) 
seal stiffness at zero eccentricity 
dimensionless stiffness 
seal lengtli 

total mass flow tlirough seal 
dimensionless mass flow 
local unit mass flow 
rotational speed 
resonant frequency 
pressure 

upstream stagnation pressure 
fluid velocity 
axial coordinate 
angle of taper 
step length ratio Lj^/L 
eccentricity ratio e/ 
angular coordinate 
dimensionless seal length fL/Cg 
fluid density 


Subscripts: 

b Bernoulli 

f friction 


L 


at z = L 
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s stepped 

t tapered 

z at location z 

1 seal entrance or large clearance region 

2 seal exit or low clearance region 

INTRODUCTION 

Annular pressure seals can significantly influence the dynamic behavior of 
rotating machinery. This is not surprising when one realizes that an annular 
clearance seal is essentially a journal bearing, although usually with larger 
clearances than normal journal bearing practice. This larger clearance re- 
duces the hydrodynamic forces generated; however, the hydrostatic pressure 
difference across the seal introduces an additional meaiis of producing seal 
forces. These forces vary directly with pressure difference across tlie seal 
(refs. 1 and 2) and can reach appreciable magnitudes in high-pressure seals 
(ref. 3). References 1 and 2 calculate stiffness and damping propei-ties of 
constant- clearance annular seals. 

Tlae variation of seal radial clearance along the axis influences the forces 
generated. The purpose of this paper is to analyze the effects of two particular 
configurations: a tapered seal in which the clearance decreases uniformly from 
inlet to outlet, and a stepped seal in which the clearance is imiform except for a 
step decrease along the seal length. Radial stiffnesses and fluid flow rates will 
be calculated using a small eccentricity analysis, and optimum configurations de- 
rived. An example of favorable influence on rotor dynamic behavior will be pre- 
sented. 

ANALYSIS 

The configuration to be an.alyzed is that of an aimular clearance seal in 
which the clearance varies along the length of the seal. The solution will be 
obtained for two clearance distributions: (1) a constant taper in the axial 
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direction; (2) two constant clearance regions wMcli join at a step. Solutions 
maj'^ also be obtained for other configurations in a straightfonvard inanner. 

The two designs considered are illustrated in figure 1. 'i?he inlet clearance h^ 
is, in general, larger than the outlet clearance h(j. Around the circumference 
the clearance is given by 

h{z, 0) = C{z) +02^ cos 0 (1) 

When the seal is eccentric, that is, when e 0, the fluid pressure in the seal 
varies in such away that a force is generated which acts to restore concen- 
tricity. To determine this force, the pressure distribution will be calculated 
and integrated over the seal area. The following assumptions are employed. 

1 . Eccentricity is small compared with the concentric clearance, that is, 

e « 1. 

2. The fluid is incompressible. 

3. Rotational effects within the seal are neglected. It can be shown from 
the results of reference 2 that shaft rotation has a minor effect in high- 
pressure seals. 

4. Fluid flow is one- dimensional in the axial direction; that is, circum- 
ferential flow is ne^ected. For high-pressure seals and small eccentricities, 
this is reasonably true. 

5. The friction factor is constant everywhere within the seal. For tur- 
bulent flow this should introduce little error, as friction factor is a weak func- 
tion of Reynolds number. 

6. Flow entrance effects are negligible. For turbulent flow, entrance 
losses are barely detectable farther than 10 hydraulic diameters from the en- 
trance (ref. 4). 

The use of the Reynolds lubrication equation, as hi a typical journal bear- 
ing, is inappropriate here because (a) the flow is usually turbulent and (b) effects 
of fluid mertia arc significant. However, a tractable solution maybe obtained 
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using principles developed for duct flow (ref. 4) . 

The axial pressure gradient due to friction is given by 


PV^ - 2f pV^ 
dzL d 2 h 2 


where the hydraulic diameter d = 2h. The mass flow per unit width (which 
does not vary with z) is 

m pVli 


Thus 


dp\ _ fm^ 
^7f ph'^ 


( 3 ) 


In addition to the friction pressure drop, die pressure will change due to changes 
in velocity (Bernoulli effect). Assuming stagnation conditions upstream of the 
seal, 




( 4 ) 


The local pressure is then given by the upstream pressure plus the pressure 
change due to velocity increase plus the friction pressui’e change. 



(5) 


For a given clearance distribution, the integral may be evaluated and the pres- 
sure distribution determined. 

For the tapered seal the clearance is given by 

h = h^ - Q!Z 

where 
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+ e cos 6 


— dz - — 


2a n. ,2 

(hj^ - a;z) hj^ 


For the stepped seal, 




Z 1 


Xj- 2 — Xj-, 

I =_i = i L-<z^L (6b) 

u3 ,,3 ^ 

hi hg 

Either of equations {6) may be substituted into equation (5) and the flow found by 


setting z = L {where h = hg) . 


,, 2 ft-0 - ly 

— + n. 


The pressure within the seal is now known, and is given by 


P=P0- (Pq- Pj) — 


— + fl-r 

<p Xj 


The radial restoring force for aii eccentric seal is found by integTating the pres- 
sure over the seal area, according to 

F=-— / / p cos 0 dz dfi (9) 

2 Wo yJo 
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The first Pq term in equation (8) does not contribute to the restoring force, 
since its integral is zero. The z integration of the remainder of equation (8) 
is straightforvvai'd. Defining 



(P - Po)^2 


it may be readily verified that 




Pq" Pl 


(o;+f)h^-fl| L’ 


(l+ljh^h2(h2-h^)+fLh2 


( 10 ) 


(lla) 


Js“-<P0" Pl^ 


h^h^L^ + fli2Lj^(L^ + 2L2) + hJhgLg + flif Lg 
"?"2 thjLj) 


(lib) 


Before the integration in 0 can be performed, the circumferential variation 
of the clearance (given by eq. (1)) must be considered. Since tlie eccentricity has 
been assumed small, aiiy terms containing the square or higher powers of e will 
be neglected. Equation (1) is now substituted into equation (11), and terms col- 
lected according to the power of e contained. Equation (11) then becomesi with 
2 

e and higher order terms omitted, 


J = 


-(PO- Pl)L 


A + Bg cos 0 
E + Gg cos 0 


(12) 


The restoring force F is given by 


P-i(Po-Pi_)LD 




A + Be cos 0 
E + Ge cos 0 


cos 0 de 


(13) 


Expressions somewhat simpler to integrate than equation (13) maybe obtained 
if the stiffness at zero eccentricity is calculated instead of the force. 
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^^27r 

K = lim (Ppi - p, ) lim / (A + Be cos 0) (E - Ge cos &) cos 0 d0 

e-vO Co£ ^2 ^ 6-^0 Uo 


•'2" 2C2E 


where, for small e, 


(14) 


E + Ge cos 8 jj2 


The constants A, B, E, and G are given by 

Aj. = H(H +A) 


B^ = H+A, 


(E - Ge cos 0) 


Ej. = h|h? +A(H +1)| 


Gt " 2At 


A^ = Hy +Ay(2 - y) +A{1 - y)^H^ + (1 - y)H' 
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E 


B^ = y (3H + 1) + 3A 1 1 + (1 - y)^ (H^ - 1)_ + {1 - y)H^ (B + 3) 


where 


Eg = + 2 a|j + (1 - y)H^ 

Gg = H^(H + 3) + 6A |y + (1 - y)H^ 

H = Cj^/Cg', A = fL/Cg; and y = L^^/L 
n 

It is well Imown that / cos” 0 d0 = 0 for an 5 '^ odd n and that 
/^27r 

^ cos^ 0 d0 = IT. Thus, equation (14) becomes 


(15) 
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^_ (Pq-Pl)LD ^ 
^2 2E^ 


(BE- AG) 


( 16 ) 


It may be noted that equation (16) agrees with the results of reference 1 when 
H = 1 and the present assumptions are applied to the analysis of reference 1. 
it is convenient to present results for a non-dimensional stiffness K where 


KCo 


K =• 


_ ir 


(Pq - Pl^^^ 2E^ 


(BE - AG) 


(17) 


Seal leakage flow may also be expressed in dimensionless terms. The total 
(dimensional) mass flow is denoted by M where 



For the concentric seal, 


(18) 


M = TrDm (19) 

The dimensionless mass flow, from equations (6), (7), (15), and (19) is 

M = 

TTDCg 

RESULTS 

The preceding expressions have been used to optimize the seal clearance 
ratio H and (for stepped seals) the step location y. The optimization was 
carried out in two ways: first to maximize the seal stiffness K, and next to 
maximize the ratio of stiffness to flow rate. The latter optimization yielded 
stiffnesses only slightly lower than maximum, while reducing the leakage flow. 
The optimization technique was similar to that employed in reference 5, and 
utilizes the Newton- Raphson method. 


M 


/(Pq - Pl^p 



( 20 ) 
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Results are presented as functions of the dimensionless seal length 
A = fL/ Cg . Figure 2 shows optimum clearance ratio for the tapered seal, and 
figure 3 for the stepped seal. For both seals, the optimum clearance ratio is 
higher for maximum K than for maximum K/M. The bto values converge at 
low values of seal length. It can easily be shown that as seal length approaches 
zero, the optimum clearance ratio approaches 2 for the tapered seal and 1. 5 
for the stepped seal. The same value of H maximizes both K and K/M for 
very short seals. 

The clearance ratio for stepped seals exhibits some unusual behavior which 
can be understood better by also considering the results for optimum step lengch 
y = Lj^/L, figure 4. The optimum step length became 1 for low seal lengths. 
This means that the large- clearance portion of tlie seal extended throughout the 
entire seal length with an .annular orifice at tlie outlet, y - 1 is, of course, the 
upper bound for step length, and when the optimum step length is 1 there is 
essentially only one parameter remaining to optimize: H. Thus, this is a dif- 
ferent regime of the solution than when y can vary freely. The cusps in the 
curves of figure 3 occur where the two solution regimes meet. 

A step length of l would not be a practical design. Thus, a second optimi- 
zation was performed in which y was arbitrarily limited to 0.8. Tliis became 
the optimum value for the entire range of seal lengths considered when maxi- 
mizing K; for maximum K/M, y == 0.8 is optimum only up to A = 0. 6. Fig- 
ure 3 shows that the optimum clearance ratio is lower when step length is re- 
stricted. 

Stiffness results . - Figures 5 and 6 illustrate the stiffnesses achieved with 
optimum tapered and stepped seals. On both figures, the stiffness of a straiglit 
seal with the same minimum clearance {calculated by setting H = 1) is also 
plotted. Figure 5, for a tapered seal, shows that the optimum seal stiffness is 
always at least 70 percent greater than that of a straight seal. Tlie stiffness 
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gain with the tapered seal is tmich greater over most of the range of seal 
lengths. Dimensionless stiffness drops as length increases. However, the drop 
is gradual enough that greater dimensional stiffness can always he obtained by 
malting the seal longer. 

Figure 6, for stepped seals, has cusps and branches in its curves corre- 
sponding to those in figures 3 and 4. A comparison of figures 5 and 6 shows 
tliat die optimum stepped seal is always stiffen than the optimum tapered seal, 
if the step lengdi is allowed to vary freely. This is to be expected, since, in 
the stepped seal, two parameters may be varied {clearance ratio and step length) 
versus only one (clearance ratio) in the tapered seal. Wlien the step length is 
restricted to a maximum of 0.3, however, the optimum tapered seal is stiffer 
for seal length A < 0.5. The stiffness of a stepped seal varies from 1.7 to 14 
times that of a similarly dimensioned straight seal, over tlie range of seal 
lengths investigated. 

Flow rates . - Flow rates through tapered seals are shown in figure 7 and 
through stepped seals in figure 8. The flow through a straight seal is shown 
for comparison in both figures. Flow through either tapered or stepped seals 
is only moderately higher than through straight seals. The relative increase 
in flow through high stiffness seals compared to straight seals rises with in- 
creasing seal length and reaches approximately 60 percent for the longest seals 
investigated. 

A comparison of figures 7 and 8 shows that, when K/ M is maximized, 
the flow rates are nearly identical for tapered and stepped seals. If only K 
is maximized, and the step length allowed to vary freely, the stepped seal 
allows a noticeably higher flow. 

Effect of variable friction factor . - The seal clearance changes appre- 
ciably from one end of the seal to the other. If this change produces a variation 
in friction factor through the seal, then results calculated assuming a constant 
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friction factor could be in error. Two flow regimes need to be considered. In 
the first, that of laniinar and transition- turbulent flow, friction factor depends 
primarily on Reynolds number. In a channel, which the seal flow passage is 
assumed to be, Reynolds number does not vary as long as the luiit mass flow 
rate, density, and viscosity are constant. Since these conditions have been 
assumed, Reynolds number, and hence friction factor, do not vary along the 
length of the seal. 

The second flow regime is that of fully turbulent, or inertial flow. For 
Reynolds numbers sufficiently high, friction factor depends only on the rela- 
tive roughness of the flow bounding surfaces (ref. 4). An expression relating 
friction factor to relative rou^iness is presented in reference 4. 

4 log, -^ + 2.28 

>fl 

This expression is rather cumbersome. However, a good approximation is 

/d\-0.29 
f = 0.036(-2-\ 

For the stepped seal, the friction factors for the two regions are then related by 

f = f ‘ ^ ^ 

The results of the ANALYSIS section are readily modified to accept different 
friction factors for the two seal sections. This was done, and the optimization 
procedure carried out. Optimum seal parameters, as well as stiffness and flow 
rate varied no more than 3 percent from the constant friction factor case. Thus 
very little accuracy is lost by assuming a constant friction factor. 
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EFFECT ON ROTOR BEHAWOR 

The rotor considered in this example is from the space shuttle main engine 
high-pressure fuel tUrbopump. Figure 9 shows the rotor. Three centrifugal 
pump stages and two axial turbine stages are all between the ball bearing sup- 
ports. The bearing span is 580 millimeters; tlie rotor weighs 56 kilograms 
(ref. 6), Annular seals are between the pump stages. The fluid pumped is 
liquid hydrogen; pressure difference across each seal varies approximately as 
the square of the rotor speed and is 136 bar at the design speed of 37 000 rpm 
(ref. 7). Seal diameter is 76 millimeters; length, 43 millimeters; and radial 
clearance, 0.13 millimeter. 

Critical speeds with original straight seals. With these dimensions and 
pressures, the flow is fully turbulent; thus the friction factor depends only 
on the relative rou^iness. For a surface roughness of 0,8 micrometer, which 
is a common machining specification, the relative rougliness d/e^, is 300 and 
the friction factor 0. 007 . Tlie dimensionless seal lengQi A is then 2. 4, From 
figure 6 for a smooth seal, K ^ 0.11. 

Critical speeds for the rotor wer-*- calculated, using a criticai speed com- 
puter program which allowed bearing stiffness (in this case, seal stiffness) to 
vary with rotor speed (the program was the same as that used in ref. 6). The 
first and second critical speeds were found to be 9250 and 17 97 0 rpm, with 
ball bearing support flexure stiffness of 30 MN/m. Rotor mode shapes are shown 
in figure 10. Since the normal speed range of this pump is from 22 000 to 
37 400 rpm, at first glance the critical speeds seem well disposed. As the 
rotor is liglitly damped, however, it is susceptible to various instabilities in 
which the exciting forces often occur at a frequency half or less that of the rotor 
speed. An example of this type of excitation is unbalanced turbine (Alford) forces 
(ref. 8) . Thus it is desirable that any resonance occur at a frequency more than 
half the rotor speed. 
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High stiffness seals . - The first critical speed is obviously the one that 
needs the largest increase. Figure 10 shows that at this critical speed, am- 
plitudes are quite high at the seal locations. Thus increased seal stiffness 
should result in an increased critical speed. To investigate the effectiveness 
of higii stiffness seals, we will assume that the existing straight seals will be 
replaced by stepped seals optimized to maximize K/M. The present overall 
seal dimensions and minimum clearance will be retained. 

The dimensionless seal length is thus the same as calculated for the straight 
seal; A =2.4. From figure 3, the optimum clearance ratio is H=1.59, From 
figure 4, the optimum step length is Y = 0.70. Figure 6 shows the stiffness K 
to be 0.29, or nearly three times the stiffness of the straight seal. With the 
stepped seals, the first two critical speeds are 10 060 and 20 190 rpra. Thus 
the first and second critical speeds have increased 9 and 12 percent, respec- 
tively. 

The second critical speed is now more than half the design speed. However, 
the first critical speed is still less than one- third of the design speed. At this 
point, it is necessary to consider the difference between a critical speed and a 
resonant frequency for a particular Operating condition. By definition, a critical 
speed occurs when the operating speed of the rotor coincides with one of its 
resonant frequencies. If the bearing and seal stiffnesses are invariant witli 
speed, and gyroscopic inertia is neglected, then the resonant frequencies are 
also invariant. For the example considered, however, the seal stiffness in- 
creases as the square of the speed (because of the increase in pump pressure 
rise); the rotor also has a large polar moment of inertia. Tlius the resonant 
frequencies will vary with rotor speed. In order to investigate the effect of 
speed on resonant frequency, the critical speed computer program was modi- 
fied to calculate resonances rather than critical speeds; the results are shown 
in figure 11. 
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The first two resonant frequencies both increase with speed. Critical 
speeds occur where the resonant frequency curves intersect the synchronous 
line (Nj^ = N) . The half- synchronous line, = N/2 is also shown; to meet 

our goal, the first resonance line should be above the half- synchronous line 
tliroughout the operating speed range. As figure 11 shows, this does not 
occur: rotor resonant frequency is less than one-half synchronous frequency 
for rotational speeds greater than 27 500 rpin. Thus, simple replacement of 
straight seals by similarly-dimensioned high- stiffness seals affords some 
improvement but does not meet our goal in this example. Additional means to 
increase the first resonant frequency include (1) changing the stiffness of tlie 
ball bearing support flexures; (2) increasing the seal length or diameter (stiff- 
ness is proportional to seal area) . However, investigation of these effects is 
beyond the scope of this paper. Stiffening the rotor is likely to have only a 
small effect, as rotor bending is very slig|lit with stiff seals at higli speeds. 

CONCLUDING HEM ARKS 

An analysis has been presented to calculate the radial stiffness of annular 
seals whose clearance varies in the flow direction. Seal parameters were op- 
timized, and results presented for seals with (1) a imiform taper in the flow 
direction; (2) two constant- clearance regions which meet at a step. Over the 
range of dimensionless seal lengtiis investigated, the optimized seal stiffness 
ranged from 1.7 to 14 times that of a similarly-dimensioed straight seal. 

High- stiffness seals can be used to raise critical and resonant speeds of rotors, 
thereby improving rotor response characteristics in some cases. 
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Figure 1. - Tapered and stepped annular seals. 
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Figure 2. - Optimum clearance ratio of tapered seat. 



Figure 3. * Optimum clearance ratio for rtepped seal. 
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Figure 10. - Critical speed mode shapes with straight seals. 
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Figure 11. - Rotor resonant frequencies with optimum stepped seals. 
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